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ABSTRACT 

Observations of type la supernovae at high redshifts have become an important tool for studying the 
geometry of the universe. The relation between the duration of the peak phase of a type la supernova's 
lightcurve and its luminosity (broader is brighter) forms the cornerstone of this measurement, yet it is 
a purely empirical relation. In this paper we show that the relation is a natural consequence of the 
radiation transport in type la supernovae and suggest constraints on the nature of the explosions which 
arise from our interpretation of the observed relation. The principle parameter underlying the relation is 
the mass of radioactive 56 Ni produced in the explosion. The relation is shown to be relatively insensitive 
to most other parameters in Chandrasekhar-mass explosions. Cosmological results are thus unlikely to 
suffer from systematic effects stemming from evolution in the explosions' progenitors. 

Subject headings: supernovae:general, cosmology:distance scale, radiative transfer 

1. INTRODUCTION 

Because they are among the brightest optical explosions in the universe, the promise of employing type la supernovae 
(SNe la) as cosmological probes has long been recognized (Kowal 1968). In the past few years, this promise has been 
realized by two groups (Schmidt et al. 1998; Perlmutter et al. 1999) who have used observations of SNe la to determine 
luminosity distances to supernovas at redshifts beyond 2 = 1. The results have been surprising: ours is a universe which 
is expanding at an accelerating rate. In several ways this result is in accord with our expectations. The observations 
imply that the age of the universe is near 14 Gy, comfortably older than the age of its contents. They also imply that the 
density of ordinary matter is well below the critical density, near £l m atter = 0.3, and this too is in agreement with most 
dynamical measurements Perlmutter et al. (1999). Perhaps the most exciting result, and certainly the most unnerving 
and controversial, is that the geometry of the universe is today dominated by an energy density which does not arise from 
ordinary matter but which is best fit by a vacuum energy density for which there are no satisfying explanations (Riess 
et al. 1998; Garnavich et al. 1998). 

To be really useful in the classic cosmological test of the redshift/luminosity-distance relation, the luminosity of a 
standard candle must have a small dispersion. The set of "typical" SNe la used for cosmology exhibits a range of almost a 
factor of three in peak brightness. This would not be sufficient to enable the current sample, which extends to a redshift of 
z ~ 1, to place useful limits on cosmological acceleration; the observed departure from a flat, matter-dominated universe 
amounts to three tenths of a magnitude (~ 30% in brightness). SNe la are useful as cosmological probes only because 
what variations are observed can, for the most part, be calibrated out, leading to a usefully small dispersion in corrected 
brightness of about 0.15 mag. 

Phillips (1993) demonstrated from a set of nearby SNe la that the rate of decline from peak in the optical lightcurve 
is correlated with peak brightness. Using the larger and more uniform data set provided by the Calan-Tololo supernova 
survey, Hamuy et al. (1996d) showed that this correlation exists in the B-, V-, and I-band lightcurves (though, interestingly, 
with different slopes than those determined by Phillips 1993 from the nearest supernovae). They found that not only the 
decline rate but also the overall width of the lightcurve centered on maximum light is correlated with peak luminosity. 
Using this same sample, Riess et al. (1995) also showed that the lightcurve shape could be used as a predictor of luminosity. 
Perhaps most suggestive from a physical point of view, using a sample of high-redshift supernovae (for which very early 
time observations are more easily obtained), Perlmutter (1997) have shown that not only is there a correlation of the 
lightcurve shape with luminosity but that, near peak, the lightcurve is the same in all supernovae up to a correlated scaling 
in time and luminosity, although quantitative errors for this relation have yet to be published. 

While this width-luminosity relation (henceforward WLR) is very well substantiated observationally, it is an empirical 
relation which has remained largely unexplained. Hoflich et al. (1996) showed that many current explosion models come 
close to this relation and suggested that the relation was determined primarily by the precise nature of the explosion 
itself. In this paper we provide the rudiments of a physical explanation for the WLR which is based upon the physics of 
radiation flow in SNe la, one in which the details of the explosion itself are largely masked. If true, this argues strongly 
against any significant bias arising in the WLR due to evolutionary effects in the supernova progenitors. 
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SNe la are widely believed to result from the thermonuclear incineration of an accreting carbon/oxygen white dwarf. 
Beyond this, however, little is known with any degree of precision. Present uncertainties include the nature and evolution 
of the progenitor system, the mass of the dwarf at ignition, and the physics of the subsequent nuclear burning. Given this 
state of our knowledge, it would at first blush appear hopeless to explain the WLR from physical principles. It may be, 
however, that the calibrating relation arises naturally from conditions present in many or all such explosions of accreting 
dwarfs and that the details of progenitor evolution and burning physics are of little consequence to the WLR, though 
fascinating in their own right. 

For want of any more certain starting point, we adopt this point of view and examine how far it may take us toward an 
understanding of the WLR for SNe la. In the first section, we reprise some of the physics of radiation transport in SNe la 
explored in two previous papers (Pinto & Eastman 2000b, a, hereafter PE I & II). We then develop a schematic description 
of the lightcurve physics which illustrates the two main properties of the explosion which determine the WLR: the mass of 
radioactive 56 Ni produced in the explosion and the rate of change of 7-ray escape. In the following sections we then present 
a sample set of lightcurve calculations which demonstrate the essential correctness of this view, and demonstrate that the 
observed WLR can be reproduced with this physics. We conclude with speculations on the probability that evolutionary 
effects can materially alter the relation observed in the local universe and lead to false cosmological conclusions. 

2. THEORY 

The WLR is the observation that the greater the luminosity at maximum light, the greater the width of the light 
curve, or, equivalently, the slower the rate of decline from maximum. As discussed in PE I, the lightcurve of a SN la 
is determined by a competition between three effects: the deposition of energy from radioactive decay, the adiabatic 
conversion of internal energy to kinetic energy of expansion, and the escape of internal energy as the observed lightcurve. 

Because the initial configuration is compact and at high density, the explosion is at first tremendously optically thick. 
The timescale associated with escape of internal energy is very long, and virtually all of the energy liberated in the 
explosion is converted to kinetic energy on a timescale approximately equal to the elapsed time. Were it not for the 
continuous energy input from decay of 56 Ni and 56 Co produced in the explosion, SNe la would be both faint and short- 
lived. Fortunately, 56 Ni decay yields an energy of 3.0 x 10 16 erg g _1 and 56 Co decay yields a further 6.4 x 10 16 erg g _1 . 

As the star expands, the optical depth decreases and the timescale over which this energy can escape declines. While 
trapped, radiant energy is constantly being converted into kinetic energy; the longer the radiation remains in the supernova, 
the less energy will escape to power the lightcurve. With a constant energy source, the luminosity would increase 
monotonically, asymptotically to balance the deposition at zero optical depth. In a real supernova, however, the deposition 
declines as radioactive material decays away and as the 7— rays, in which the decay energy is emitted, escape. Peak 
luminosity thus occurs roughly when the product of the rising escape fraction and the decreasing deposition reaches a 
maximum. This is complicated somewhat by the fact that the balance is not instantaneous; the energy emerging as UV 
to IR photons at any time near peak light has been deposited over a significant fraction of the elapsed time. Arnett 
(1982) first demonstrated that at maximum light the instantaneous bolometric luminosity is approximately equal to the 
instantaneous rate of energy deposition by radiative decay. Earlier, the expansion timescale is short, the escape time is 
long, and the deposition rate is high. More internal energy goes into expansion than escapes, and a store of trapped energy 
is built up. Subsequent to maximum light, more energy escapes than is converted to kinetic energy, and the surviving 
store of trapped energy is released faster than the deposition. 

The WLR is thus seen to be a relation between the maximum-light luminosity and the radiation escape time; the 
brighter the supernova, the longer its escape time. Radiation escapes from an optically thick medium by diffusion, and 
we can tentatively identify the escape time with the diffusion time - the average time it takes a photon to random walk 
its way to the surface, given by the number of scatterings times the time of flight between scatterings. 

Following the usual random walk argument, the diffusion time can be written as 

t d = ^ (1) 

where R is a characteristic size of the system (in this case, the radius of the ejecta), r e ff is the effective optical depth, 
and a is a constant of order unity. By the argument above, if the expansion (elapsed) time and escape (diffusion) times 
are about equal at peak, we can write 

aRT e ff CtVtTeff 

t d = t= — = — (2) 



which gives a diffusion optical depth at peak of 



10 9 cm s- 1 



Teff « 30 ^ j (3) 

(taking a ~ 1). This value is far smaller than the monochromatic optical depth to most trapped photons, which is typically 
nearer to 3000 for peak conditions in a SN la. The resolution to this seeming paradox, of course, is that the typical photons 
in an atmosphere do not carry the flux - radiation escapes where it can, at energies where the monochromatic optical 
depth is smaller. 

At stellar densities, there is tight collisional coupling between the radiation field and the plasma, keeping the local 
photon distribution function at its blackbody limit. While the bulk of the energy density exists as UV photons, these 
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photons are trapped by the high UV opacity. The bulk of the flux is transported at longer wavelengths on the Rayleigh- 
Jeans tail of the distribution, where the optical depth is typically less, "side-stepping" the higher UV opacity. The lower 
density of these photons is more than overcome by the decreased diffusion time they suffer, and the appropriate opacity 
to use in describing the flux is the familiar Rosseland mean. 

For this picture to work, the relaxation rate of the trapped photon distribution toward a blackbody must be short 
compared with the rate at which energy escapes from the long-wavelength tail of the distribution. This relaxation is 
usually mediated by collisions, and its rate must therefore scale as the square of the density. If the optical depth to the 
majority of photons remains large (so that radiation remains trapped) while the density decreases, eventually collisions 
will not be able to maintain sufficient energy density in the long-wavelength tail of the photon distribution where the 
optical depth is less; the system will depart from local thermodynamic equilibrium (LTE). 

PE II demonstrated that this has indeed become the case in SNe la by maximum light. The escape time is still mediated 
by the rate for converting energy from short-wavelength, UV photons into long wavelength, OIR photons. The difference 
from the more typical LTE situation arises from the nature of opacity in rapidly-expanding media and the character of 
atomic physics in iron-group ions. 

In the homologous expansion which develops soon after explosion, the velocity gradient is isotropic. Much as in the 
expanding universe, there is a one-on-one mapping between path-length and redshift, and a photon traversing the ejecta 
will accumulate a red shift in the local gas rest frame which will bring it into resonance with line transitions. One can 
show (Eastman & Pinto 1993) that under the conditions in SNe la, the time a photon spends scattering in a line is very 
short compared with the time it spends moving between lines. For the purpose of estimating the escape time, then, each 
line traversed thus acts as a single scattering, and the distance measured in units of a photon mean free path (the optical 
depth) will be just the number of optically thick lines traversed in accumulating a given Doppler shift; the corresponding 
effective opacity is known as an "expansion opacity" (cf. Eastman & Pinto 1993; Wagoner et al. 1991; Karp et al. 1977; 
Castor et al. 1975). 

Using simple analytic estimates, PE I showed that the central temperature in a M c h explosion at t < 20 days is 
T > 13, 000 K and the peak of the blackbody spectrum is in the UV near A < 2200 A. For an iron-peak composition in 
the relevant range of temperature and density, the spectral density of optically thick lines in the UV is very large and 
is a strongly decreasing function of wavelength from the UV through the optical (cf. Figure 2 of PE II). As long as the 
Sobolev optical depth in these lines remains large, the diffusion optical depth will not change as the density decreases, 
allowing a departure from LTE not seen at depth in static atmospheres where the optical depth is linearly dependent 
upon density. The large spectral density of lines at UV wavelengths is far too great for there to be significant transport at 
UV energies; the time it takes a UV photon to random walk to the surface is long compared to the expansion time. This 
is substantiated by the observation of very little UV flux in maximum light SNe la in spite of their high temperatures. 

If the Rayleigh-Jeans tail of the photon distribution (at wavelengths where significant transport can take place) could 
be re-populated by a mechanism more efficient at low densities than collisional excitation, significant flux could still be 
developed, leading to an escape time characterized by the required effective optical depth near 30. 

When a photon is absorbed in a line, the possible outcomes of this interaction include: 1) re-emission in the same 
transition, followed possibly by re-absorptions and emissions - a cycle repeated many times until escape from the red 
wing of the line; 2) collisional destruction, transferring the photon energy by Coulomb interactions to a free electron; and 
3) fluorescence, where the excited state populated by absorption de-populates by radiative decay to a state different from 
the original. In general, the electron density is too low in maximum light SNe la for collisional destruction to have much 
effect. Thus, a UV photon absorbed on the high energy (blue) side of a line will either be re-emitted on the red side, or 
fluorescence will occur, with the photon replaced by several longer- wavelength photons (usually), but occasionally by a 
shorter wavelength photon (less often). The efficiency for this conversion process obviously depends on the number and 
strength of decay channels which would result in the emission of an OIR photon. For the UV radiation trapped in SNe la, 
fluorescence is the most likely outcome. Fluorescence is much more important in SNe la than in other astrophysical 
objects because they are composed in large part of iron-group elements with enormously complex atomic structures and 
very high spectral densities of line transitions in which fluorescence can occur. 

A key point is this: the number of OIR transitions is greatest in the neutral and low ionization states, and decreases 
with increasing ionization. Thus the UV^OIR fluorescence process is more efficient in Co II than in Co III, and more 
efficient in Co III than in Co IV. Consequently, the conversion process is sensitive to the ionization. In the approximate 
calculations reported below, we assume LTE ionization and excitation despite densities which are too low for this to be 
really true. We shall thus use temperature as a proxy for ionization, but the reader should keep in mind that this is an 
approximation, and it is actually ionization which is the controlling variable. 

In light of the importance of fluorescence, we can sketch the escape of radiation from SNe la in the following fashion. 
Imagine a store of deposited energy, trapped in the form of UV photons. These UV photons are created by the absorption 
and thermalization of decay 7-rays and removed by conversion to kinetic energy and to OIR photons. For the energy 
density of UV photons, Ejjv , we can write the following schematic rate equation: 

dE uv ■ E uv 

—gf- = pSde P (t) - -j- - r f E uv (4) 

where pSde P (t) is the time-dependent deposition rate, Ejjv/te is the rate of conversion of radiation to kinetic energy, and 
rfEuv is the rate for converting UV to optical and infrared energy. We have assumed here that the diffusion time suffered 
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by the UV energy is infinitely long; transport in the UV is unimportant. Similarly, for the OIR photons one can write 

dEom E ir E ir „ , . 

-8T = — — + rfEuv (5) 

where td is the typical OIR diffusion time. Here we have ignored the possibility that deposition results in the direct 
emission of OIR photons. If optical photons are processed on a short enough time scale - for instance, if t d < t e - then 
we can set dEom/dt ~ and solve for Eoir- 

Eoir = —7- — —fy. — : Euv- (6) 

l/t e + l/t d + r f 

Finally, by adding equations (4) and (5) together, and making use of equation (6), we obtain an equation for the total 
energy density, E tot : 

dE to t a / . \ E t ot Etot /„n 
~dT = pSd ^ ~ — ~ 17 (?) 

where the effective diffusion time, td, is defined as 

- _ t d + t e + rftdt e . . 

td = — ^ — • (8) 

Written this way, it is easy to see that as rf - the UV to OIR fluorescence rate - decreases, the effective diffusion time 
goes up. In terms of the opacity r/ can be written in LTE as 

where B v is the Planck function at the local gas temperature, and voir is a frequency which approximately divides the 
UV from the optical. The first fraction is just the Planck mean opacity in the UV, and the second fraction is the fraction 
of absorbed energy which is radiated at optical wavelengths. Since, for instance, Co II has a larger opacity in the optical 
than Co III, increasing the ionization from once- to twice-ionized will decrease the fraction of absorbed energy which is 
re-radiated in the optical, making r/ smaller, and therefore increasing the effective diffusion time. 

We can now understand the underlying physics of the WLR: as the amount of 56 Ni is increased, the interior temperature 
and ionization go up, which decreases r j and increases td- Stated another way, brighter SNe la (i.e. ones with more 56 Ni) 
are hotter inside, and therefore have longer diffusion times (i.e. evolve more slowly through peak). 

This cannot be the whole story, however. As long as most of the 56 Ni is unmixed and its mass fraction X 56 ss 1 in the 
56 Ni-dominated layers, increasing the total amount of 56 Ni in the star does not increase the local heating rate per gram, 
and therefore would not, by itself, result in a temperature increase. Overlooked so far is the fact that, if the additional 
56 Ni is added at higher velocity, then there will be an increase in the effective optical depth. That is, in stars with more 
56 Ni (in a layered composition), the velocity of the outer edge of the 56 Ni is moving faster than is the outer edge of the 
56 Ni in a star with less 56 Ni. The higher the velocity of the outer edge of the 56 Ni, the more lines through which a photon 
must Dopplcr shift before escaping. This represents an optical depth increase which would increase the escape time and 
therefore make the gas hotter. Putting iron-peak elements at higher velocity is, by itself, not sufficient to increase r e // 
enough to fully explain the observed correlation between light curve width and luminosity. Rather, the slight increase in 
optical depth provided by putting 56 Ni at higher velocities is amplified by its effect on temperature and ionization, and 
thus the escape time. It is difficult to cleanly separate these two effects, but in the next section we present results which 
strongly suggest that ionization is the dominant mechanism determining the escape time. 

It is worth stressing here that fluorescence is not a one-way process. Energy can be transported to shorter wavelengths 
by "combining" photons. A photon can be absorbed, raising an ion to an excited state. The ion can then absorb another 
photon into this state, raising it to yet a higher-energy state from which a UV photon can be emitted. From obvious 
thermodynamic considerations (embodied in Rosscland's Theorem of Cycles) this "combining" of photons must in general 
be less frequent than "splitting" of UV photons. Near the peak in the lightcurve, however, over much of the volume of 
the supernova these processes are nearly in equilibrium; it is the net rate of transport of energy to longer wavelengths 
which determines the escape time, not merely a one-way fluorescent cascade. Indeed, Mazzali (2000) has suggested that 
"combining", in lower-optical depth regions above the photosphere, is a significant contributor to the UV emission at 
maximum light. 

Another avenue for time-dependence to enter the equations above is through the energy deposition term Sd ep (t)- The 
deposition is determined both by how much 56 Ni is present and by what fraction of its decay energy is deposited in the 
gas. Because the Compton opacity to MeV 7-rays is mainly absorptive, the deposition fraction is roughly proportional 
to 1 — exp(—Tj), where t 7 is the mean gamma-ray optical depth and is in turn roughly proportional to the mean column 
depth to 56 Ni. More-rapid expansion will lead to a lower mean optical depth to 7-rays at a given time, as will mixing 
radioactive material to higher velocities. Lower column depths will in turn lead to a more-rapid decline in the deposition 
fraction. If S dep (t) declines more rapidly, the lightcurve will peak earlier. An earlier peak will have more 56 Ni still present, 
and this will tend to offset the larger escape fraction. The lightcurve will thus evolve more rapidly, leading to a narrower 
lightcurve for it luminosity. 

We have thus identified three properties of the explosion which may have a strong effect upon the WLR. The mass of 
56 Ni determines the ionization and through it the mean escape time. The expansion velocity (or specific kinetic energy) 
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determines the rate of decline of the energy deposition. The radial distribution of 56 Ni can have an effect upon both the 
escape time of the thermalized radiation and the time evolution of the deposition. In the following sections we present a 
small set of synthetic lightcurves which illustrate these effects and their relative importance. 

3. COMPUTATIONS 

A complete model for supernova lightcurves would begin with a model for the evolution of the progenitor. This would 
be evolved to ignition, the subsequent explosive burning followed to completion, and a transport calculation performed to 
determine the lightcurve. Unfortunately, a fully predictive model is not yet available from such a procedure. While there 
are schematic models for all of these stages, it is not yet possible to simulate SN la explosions without introducing very 
simple parameterizations to stand in for very complex physics we do not yet understand. 

The evolution of the progenitor is undoubtedly the most uncertain part of this process, yet there is the possibility that 
"convergent evolution" of some sort, perhaps ignition near the Chandrasekhar mass in single degenerate systems, reduces 
the sensitivity to initial conditions. On the other hand, the initial conditions in the various double-degenerate scenarios 
might well be the single most important piece of physics determining the subsequent explosion. In the burning phase, the 
physics of turbulent combustion and the possible spontaneous transition to detonation are probably the most important 
and least tractable effects. While progress is being made on all of these fronts (Reinecke et al. 1999; Iwamoto et al. 1999), 
it will be a long time before truly predictive models of progenitor evolution and explosion become available. 

Our aim in this work is to explore the radiation transport effects which influence the WLR. We have thus chosen to 
construct a very simple model for the structure of the explosion. We start with model DD4 of Woosley & Weaver (1991), 
an M c h explosion which burned 1.26Mq to Si or heavier elements to yield 1.2 x 10 51 ergs of kinetic energy, producing 0.63 
M Q of 56 Ni. We have modified the mass and distribution of 56 Ni in this model to explore the effects such modifications 
have on the lightcurve. The density as a function of velocity in the original model DD4 is preserved in all models. The 
structure of DD4 is presented by Woosley & Weaver (1991) and by PE II. 

An initial grid of models with different 56 Ni masses was obtained by mapping between 56 Ni and 28 Si according to 

X 56 (M) = X 5 6 (M) exp (-0(M)) 

X 2S (M) = X 2 ° 8 (M) + X° 56 (M), (1 - exp (-#M))) 

where 

0(M)=max(a(M-M cot ),O). (10) 
To increase the 56 Ni mass, all isotopes below M — M cut were mapped to Ni: 

X m (M) = 1 - H(M - M cut ) + X° 6 (M)H(M - M cut ), (11) 
where H{x) is the Heaviside function, and for A ^ 56, 

X A {M) = X° A (M)H(M - Met). (12) 

This simple scheme is actually well-motivated in a physical sense. The total energy released from burning a C/O 
mixture to intermediate elements in the Si group is very nearly the same as for burning to iron-peak elements. The 5 Ni 
mass thus need have little necessary connection to the explosion kinetic energy. The most stringent constraints on the 
distribution of elements in velocity come from peak-phase spectroscopy, though this analysis is hampered by our inability 
to compute time-dependent NLTE models. 

As discussed in PE I, the initial thermal energy in the explosion goes entirely into kinetic energy via the second term 
on the r.h.s. of equation (4) and hence has no effect upon the observable part of the lightcurve. All the models presented 
here were started with the same initial internal energy. 

The lightcurves were computed as described in PE II, using time-dependent, multi-group radiation transport with an 
LTE equation of state and including the elements He, C, O, Si, S, Ca, Fe, Co, and Ni in ionizations up to Ni XIII. 6000 
energy groups were employed, providing moderate resolution; doubling the number of groups does not significantly alter 
the results, though halving it does leads to marked effects. We are forced to employ LTE level populations because of 
computer-time constraints - a NLTE calculation would take as long per time step as a whole lightcurve in LTE (several 
days on a fast workstation). All current NLTE calculations (Hoflich 1995; Nugent et al. 1997; Mazzali & Lucy 1993) 
employ temperature structures which are not consistent with the history of the ejecta's evolution and therefore cannot 
be expected to give correct time-dependent spectra and lightcurves, as discussed in PE II. While these LTE calculations 
obviously suffer as well, the ability to include time-dependent physics correctly is evidently necessary to study the effects 
discussed above. 

The WLR is an observed relation between peak magnitude and a measure of width. To obtain an accurate result in 
the face of measurement error and the temporal sparsity of observed data, both the peak magnitude and the width are 
measured by fitting observed data to templates determined from well-observed supernovae (Hamuy et al. 1996b; Riess 
et al. 1995). It is important to employ the same methods to measure the peak magnitude and width in our models. 
We "observe" our synthetic lightcurves by employing the MCLS (Multi-Color Lightcurve Shape) method of Riess et al. 
(1995), employing the revised vectors used in Riess et al. (1998) (kindly provided by Dr. Riess). We have also determined 
the Amis(-B) decline-rate parameter for our models using the x 2 template-fitting technique described by Hamuy et al. 
(1996c). The results of fitting from both of these techniques are shown in table 1. 

We note that the observational difficulty of obtaining early lightcurve data has hindered the construction of templates 
for times earlier than ten days before peak. Thus, by construction, template-fitting techniques are most sensitive to 
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post-maximum decline rates; their ability to discriminate between different pre-maximum lightcurvc shapes is seriously 
hampered. 

Observationally, WLR relations with similar scatter exist in the B, V, R, and I bands - color evolution during the first 
30 days is well correlated with luminosity and width. At later times, all SNe la show a remarkable homogeneity in color. 
In the absence of reddening, determining the width of the lightcurve in a single band is sufficient to uniquely determine 
a supernova's peak magnitude. The importance of fitting a supernova's color evolution is primarily in determining 
extinction. As our models do not suffer from this observational uncertainty, it is sufficient to fit to a single band. 

We remind the reader as well that individual supernovae show significant departures from the best-fit templates (cf. 
Hamuy et al. 1996a); SNe la do not form a strictly one-parameter family. While the width of lightcurves measured by the 
templates is the largest variation from object to object, there are other features in the lightcurve which are not correlated 
with the decline rate. A good example of this is the pair SNI99IT and SNI992bc which have very similar decline rates 
over the first 15 days past peak but very different rates over the subsequent 45 days (Hamuy ct al. 1996c). 

The details of the atomic physics and NLTE effects determine the color evolution of the supernova. We have assumed 
in this study that such effects are less important to the bolometric luminosity, merely reapportioning flux between the B, 
V, and R bands which comprise most of the luminosity (Suntzeff 1995). We have thus chosen to take fits to the WLR 
in the V band only as most representative of the overall relation between width and luminosity. Below we present fits to 
B, V, and R simultaneously as well, but as an accurate determination of the color evolution depends more upon specific 
details of the atomic models and NLTE effect, we regard these fits as less reliable indicators of the underlying physics 
than the trends in a single band. 

As the ejecta evolve, the effects of trapping and fluorescence become less important to the bolometric lightcurve. The 
colors at these later times, however, are very strongly influenced by NLTE effects. LTE calculations cannot be expected 
to provide an adequate representation of the lightcurve physics as the supernova becomes a nebula. Thus, we fit only to 
the peak phase of the lightcurve and estimate the sensitivity of fitting to various time intervals by fitting from 10 days 
before peak (the earliest time in the templates) until 20 and 40 days past peak. 

4. RESULTS 

Figure (1) compares the MCLS A parameter (a measure of width) to the peak V absolute magnitude determined from 
the models described above. The MCLS method parameterizes the shape of a supernova's lightcurve as 

M(t) = A(t) + A x L(t) + A 2 x Q(t) (13) 
where A(t) is, roughly, the mean lightcurve shape of a sample of 27 well-observed supernovae (Riess et al. 1998) and 
L(t) and Q(t) are functions describing a quadratic fit to the departure of individual supernovae from this mean shape as 



Table 1 
MCLS Lightcurve Fits 



model 


M( 56 Ni) 


M peak 


Am 15 a 


M peak h 


A 20 c 


M peak d 


A 40 o 
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A 20 g 


A 40 h 




(M ) 


(actual) 




(Amis) 




(A 20 ) 




(A 40 ) 


(BVR) 


(BVR) 


DD4/27 


0.27 


-18.90 


1.42 


-18.87 


0.470 


-18.98 


0.310 


-18.87 


0.509 


0.479 


DD4/38 


0.38 


-19.17 


1.47 


-19.15 


0.348 


-19.10 


0.241 


-19.21 


0.460 


0.427 


DD4/48 


0.48 


-19.34 


1.37 


-19.34 


0.121 


-19.33 


0.159 


-19.23 


0.220 


0.235 


DD4/55 


0.55 


-19.42 


1.24 


-19.43 


-0.054 


-19.51 


0.123 


-19.47 


0.160 


0.157 


DD4 


0.64 


-19.52 


1.22 


-19.52 


-0.123 


-19.58 


0.131 


-19.32 


0.160 


0.157 


DD4/90 


0.90 


-19.78 


0.87 


-19.81 


-0.443 


-19.91 


-0.542 


-20.02 


-0.704 


-0.710 


Ml 


0.64 


-19.56 


1.24 


-19.55 


-0.126 


-19.58 


0.018 


-19.43 


0.219 


0.244 


NiFe 


0.32 


-19.01 


1.75 


-19.01 


0.633 


-18.81 


0.502 


-18.94 


0.655 


0.543 


DD3 


0.93 


-19.90 


0.87 


-19.90 


0.083 


-19.36 


-0.161 


-19.62 


-0.698 


-0.525 



a Ami5 measured from the V lightcurve. 
b Peak V magnitude as measured by the A77115 fit. 
C MCLS width measured from the V lightcurve from (-10,20) days. 
d Peak V magnitude as measured by the A 2 o fit. 
C MCLS width measured from the V lightcurve from (-10,40) days. 
Peak V magnitude as measured by the A40 fit. 

g MCLS width measured from the B, V, and R lightcurves from (-10,20) days. 
h MCLS width measured from the B, V, and R lightcurves from (-10,40) days. 
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Fig. 1. — Peak V magnitudes vs the MCLS A parameter (a measure of lightcurve width) fit to V lightcurves over the interval (-f0,20) 
days (top panel) and (-f0,40) days (bottom panel) from peak. The squares correspond to calculations made by varying the 56 Ni mass in 
model DD4. The stars show the three variations from these models as discussed in the text. The diagonal curve in each panel is the observed 
quadratic relation between M^ a and A using the Cepheid calibration of SNe la. The approximate magnitude of the observed dispersion 
about this relation, corresponding to ±0.15 mag and ±0.125 in A, is indicated in the top panel 



parameterized by A which is thus a measure of lightcurve width. The quadratic relation between A and V pea k is indicated 
by the diagonal curves in the figure. The absolute luminosity of the relation is calibrated by the Cepheid distance scale 
of (Gibson et al. 2000) corresponding to a Hubble constant of 68 km s _1 Mpc -1 . The square symbols correspond to the 
models with 0.27, 0.38, 0.48, 0.55, 0.63, and 0.90 M of 56 Ni in the structure of model DD4. The star symbols represent 
other models to be discussed below. The top panel shows fits over the interval (-10,20) days relative to peak and the lower 
panel over the interval (-10,40). The quality of the fit to the mean WLR relation is somewhat better for the scaled- 56 Ni 
models for the shorter interval and is probably a consequence of the increasingly inappropriate assumption of LTE at 
later times. 

In both cases, the fits show a dispersion comparable to or smaller than the observations (~ 0.15 mag) and (remarkably 
and probably fortuitously) virtually zero departure from the absolute calibration. It is clear from this figure that the V 
WLR is fit remarkably well by the simple prescription of changing only the 56 Ni mass. 

Figure (2) shows the same lightcurves fit by the template-fitting (A 17115(B)) method of Hamuy ct al. (1996b) on the 
interval (-5,26) days. In this method, the models are fit to each of six template lightcurves spanning the range of well- 
observed supernovae. The templates are characterized by their Ami5(B) parameters (roughly the decline in B magnitude 
over 15 days from peak), while the actual epochs fit are determined by the whole peak phase of the lightcurves. The best- 
fit A 17115(B) is then determined from the minimum of a quadratic fit to the reduced x 2 as a function of the parameter. 
Rather than extrapolating, models broader than the broadest template are assigned the minimum value of A 17115(B), 
accounting for the slight pile-up of points at the bright end of the relation. 

The filled symbols are the same as those of the previous figure, and the open symbols are fits to 59 well-observed 
supernovae from Phillips et al. (1999). The curve is the mean Auii5(B)-Mv relation, on the same calibration, from 
(Phillips et al. 1999). Once again, the models fit the relation well within the scatter of observed supernovae. The 
significant difference in curvature between the mean A TOi 5 (£>)-luminosity and A- luminosity relations is notable and 
appears to be the result of the differing ways in which the lightcurve shapes are parameterized; both techniques achieve 
similar reductions in the scatter about the Hubble relation. 

Figure (3) shows the V lightcurve family from the six scaled DD4 models compared with the best-fitting MCLS 
templates. The shapes of the lightcurves depart less than 0.05 magnitude from the templates throughout the first 45 days 
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Fig. 2. — Peak V magnitudes vs the Amis(B) parameter (another measure of lightcurve width) measured from V lightcurves and on 
the same luminosity calibration as the previous figure. Open symbols are data for 59 supernovae from Phillips et al. (1999), on the same 
calibration as in the previous figure. The stars show the three variations from the scaled-DD4 models as discussed in the text. 

(about 25 days past V peak) with the exception of the 0.90 M©model, where the inflection at 32 days is not seen in the 
templates. This error is well within the dispersion of data from individual supernovae about the template lightcurves (c/. 
Riess et al. 1998). At later times, the LTE approximation becomes increasingly poor and may account for the too-shallow 
decline rates, though real supernovae show similar departures from the templates at this time. 

Figures (4) and (5) show the B, V, and R lightcurves of the 0.55 M Q and 0.90 M models, respectively. The DD4/55 
model displays a generic problem with most of the models; while the B lightcurve is the correct width, near peak it is 
too flat compared to the templates. The departure from the templates before peak is less disturbing, and is probably 
more a problem with the templates than the models, as there is comparatively little data so far before maximum. The 
IR secondary maximum is clearly seen in the R-band curves. 

We now turn to variations on these scaled-DD4 models to investigate the sensitivity of the WLR to other changes in 
the explosion. The models presented so far have nearly pure 56 Ni from the center out to a mass coordinate nearly equal 
to the Ni mass - radioactive material is concentrated as much as possible toward the center. This is a natural result of 
the explosion physics. White dwarfs are densest at the center. When matter burns explosively, the higher the density, 
the larger the ratio of the hydrodynamic timescale to the the burning timescale, and the more the distribution of burning 
products is skewed toward heavier elements. In the absence of significant hydrodynamic mixing after material has cooled 
below the burning temperature, a natural consequence of explosive burning is that the gradient of atomic weight will 
point inward, with the heaviest burning product, 56 Ni, concentrated toward the core. 

To explore the effect hydrodynamic mixing might have on the lightcurve, we have altered the distribution of 56 Ni in 
model DD4 without changing its mass. A significant mass fraction of 56 Ni extends approximately to O.85M in the original 
model (primarily as a consequence of producing O.16M of non-radioactive iron-group material). The starred points at 
My — 19.56 in Figures (1) and (2) show the result of homogenizing the composition of this model over the mass range 
0.7 to 1.1 M in model Ml. This decreases the peak bolometric magnitude by 0.1 mag and leads to significantly narrower 
B and brighter I lightcurves. The resulting V and R magnitudes are virtually identical, 0.04 and 0.02 mag brighter, 
respectively. The 56 Ni nearer the surface cooled more rapidly. The lower ionization gas converts energy from shorter- to 
longer- wavelength photons more efficiently, circumventing the effect of the higher expansion optical depth from 56 Ni at 
higher velocities. Altering the distribution of 56 Ni would appear to affect the color of the emergent flux more than its 
overall escape properties. An important conclusion is that, within the limits of fitting O.63M of 56 Ni in an M c h white 
dwarf, the WLR is insensitive to the distribution of radioactivity in the explosion. 
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Fig. 3. — The six V scaled-DD4 lightcurvcs. The solid lines are the synthetic lightcurves; the dotted lines are the best-fit MCLS template 
lightcurves fit over the interval (-11,20) days with respect to peak. 



The scaled-DD4 models traded off production of 56 Ni with intermediate mass elements. This is the only possibility 
afforded if the supernovae are to satisfy nucleosynthetic constraints. For example, the production of large quantities of 
non-radioactive iron-group elements in the majority of supernovae is firmly ruled out. Nonetheless, in model NiFe we have 
converted, in place, half of the mass of 56 Ni to 56 Fe, thus changing the heating rate per gram of the iron-rich material 
without changing its distribution, cooling, or 7-ray escape properties. The lower temperature in the 56 Ni-rich material 
leads to lower ionization and hence to more efficient escape. In spite of having 0.32 instead of O.27M0 of 56 Ni, the model 
has a much narrower lightcurve than model DD4/27, as shown by the stars near My = —19.01 in the Figures 1 and 2. 
In model Ml 56 Ni was mixed outward with overlying Si-group material. However the amount of material involved in the 
mixing was a small fraction of the total mass of 56 Ni. Consequently the effectiveness of mixing in lowering the heating 
rate per gram and narrowing the lightcurve was small. 

Interchanging Fe- with Si-group elements makes relatively little difference in the amount of energy released by burning 
a C/O mixture; at most 6.2 x 10 17 erg g _1 is released burning to Si & S versus 7.9 x 10 17 erg g _1 for burning to NSE. 
It is thus quite possible to imagine a series of explosions in which the Si/Fe ratio varies more strongly than the specific 
kinetic energy. On the other hand, in most current models a greater production of 56 Ni is accompanied by a greater total 
amount of burning, leading to higher velocities. Model DD3 of Woosley & Weaver (1991) produced 0.96 M Q of 56 Ni. It 
burned 10% more mass to the Si-group or above in total and has 10% greater kinetic energy than model DD4. In spite of 
its high 56 Ni mass, the lightcurve from this model is considerably narrower than that of model DD4/90 and is shown by 
the stars in Figures (1) and (2) at My = —19.90. The higher kinetic energy in this explosion resulted in a lower column 
depth, allowing decay 7-rays to escape more readily after peak. By 20 days, ~ 6% of the decay energy in model DD4/90 
escaped directly as 7-rays, increasing to 35% by day 40. In model DD3, the escape fraction was more than double at 
20 days, at ^14%, rising to 46% by day 40. A shortened time to 7-ray transparency leads to a shorter rise-time (by 3.5 
days), a narrower peak, and a greater decline rate in the lightcurve. The increased 7-ray escape was, however, more than 
compensated for by the shorter rise-time, and the peak magnitude is even brighter than the difference one might predict 
based solely on consideration of the difference in total 56 Ni mass. It is interesting to note that the A 17115(B) fitting 
procedure did not pick up this variation in lightcurve shape as it weights more heavily the decline rate after peak than 
the MCLS templates; the two models had very similar post-maximum declines. 
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Fig. 4. — B, V, and R lightcurves of model DD4/55 and the best-fit MCLS templates on the interval (-10,20) days relative to peak. 

While far from a thorough investigation of all the variations one can imagine on explosions of M c h white dwarfs, the 
results presented above are suggestive. They show that it is possible to create a series of models which exhibit, to the 
best of our current calculation ability yet well within observational errors, the observed WLR. 

It seems clear from the arguments of §2 and the results of these simple experiments that the 56 Ni yield and the 7-ray 
escape fraction play the determining roles in the relation between the lightcurve's width and peak luminosity. The models 
show a remarkable insensitivity to the distribution of radioactivity and to the composition of non-radioactive material, 
within limits set by achieving the appropriate explosion energies and acceptable nucleosynthesis. 

The mass of 56 Ni sets the scale of the peak luminosity and the heating rate of the gas. The resulting temperature 
determines the ionization and hence the cooling rate of the gas. More 56 Ni leads to more heating, higher temperatures, 
less efficient cooling, and hence broader and brighter lightcurves. 

The efficiency of 7-ray escape sets the evolution rate of the energy deposition. Models with higher overall velocities 
and/or with 56 Ni distributions extending to higher velocities exhibit larger 7-ray escape fractions at a given time. This 
increases the evolution rate of the lightcurve, with greater escape leading to narrower lightcurves. Models with higher 56 Ni 
yields tend to have their radioactivity distributed over a broader range in velocity, leading to larger escape fractions. This 
greater escape fraction acts to oppose the tendency for broader lightcurves resulting from higher temperatures, decreasing 
the slope of the WLR. 

We have reproduced the observed relation in V by varying only the 56 Ni yield. One solution to reproducing the observed 
WLR is thus to produce an evolution scenario which leads to supernovae with similar explosion energies yet differing 56 Ni 
yields. Because higher 56 Ni yields and higher velocities tend to produce opposing effects on the lightcurve, this solution is 
almost certainly not unique. A weak dependence of energy on 56 Ni yield will likely lead to acceptable lightcurves as well. 

While to date the models which best fit lightcurves and spectra are all at the Chandrasekhar mass (with the possible 
exception of models for SN1991bg), the masses of all SN la explosions may not be identical. Varying the explosion mass 
primarily affects the 7-ray optical depth. As long as a given 56 Ni mass is accompanied by a column depth similar to 
those exhibited by the M c h models considered here, a variable-mass explosion scenario would not appreciably alter our 
conclusions. For sub-M c h model such as those of (Livne & Arnett 1995; Woosley & Weaver 1994), the 56 Ni produced by 
detonation in surface layers does not contribute appreciably to the optical lightcurve due to the low 7-ray optical depth. 
In such models, only the 56 Ni produced by carbon detonation in the core C/O dwarf would affect the WLR. 

The exact values of the slope and calibration of the WLR obtained in this work must be, to some extent, fortuitous. 
EP II showed that the conditions experienced by matter in SNe la are not consistent in detail with the assumption of 
LTE. The shape and colors of the lightcurve depend sensitively upon the average branching ratio from UV to OIR photons 
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Fig. 5. — B, V, and R lightcurves of model DD4/90 and the best-fit MCLS templates on the interval (-10,20) days relative to peak. 

through a broad range in energies. While many strong lines are unlikely to be missing from our atomic models, all such 
models are seriously incomplete and the statistical properties of the line distributions and branching ratios are unlikely 
to be correct in detail. The basic physics, however, will not change significantly, and thus we have some confidence in the 
general trends the models indicate. 

By fitting our models using the same procedures used to calibrate observed SNe la, we have demonstrated that the 
calibration is insensitive to fairly large changes in the underlying explosions. This suggests that the small differences in 
initial conditions which might arise from evolutionary effects between the local universe and that at z ~ 1 arc unlikely 
to appreciably affect the supernova cosmology results. For example, mixing in model Ml leads to much less than 0.1 
magnitude difference from the corresponding unmixed model. Even the unphysically large change in abundances in model 
NiFe, strongly ruled out by chemical evolution considerations, leads to less than 0.2 magnitude departures from the MCLS 
and barely 0.3 magnitude departures from the A 17115(B) relations; allowable changes in relative iron-peak abundances 
would lead to insignificant effects. 

Altering the specific kinetic energy, as in Model DD3, is the only way we see to achieve significant departures from 
the WLR. Increasing the kinetic energy by 10% was seen to lead to a much narrower lightcurve (at least when fit by the 
MCLS templates), primarily by more than doubling the 7-ray escape from 6% at peak in DD4/90 to 14% in DD3. This 
change is in the wrong direction, however. To account for the 0.3 magnitude deficit in brightness seen at large redshifts, a 
shift to lower explosion energy would be required. Halving the 7-ray escape from 6% at peak to 3% is a smaller fractional 
change in deposition and would lead to a smaller change in width. A future paper will explore these effects in greater 
detail. 

We wish to thank M. Hamuy for many useful discussions and for performing the Ami5(B) fits to the lightcurves, and 
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Science Foundation (CAREER grant AST9501634), by the National Aeronautics and Space Administration (grant NAG 5- 
2798). This work was also performed under the auspices of the U.S. Department of Energy by University of California 
Lawrence Livermore National Laboratory under contract No. W-7405-Eng-48. Philip Pinto gratefully acknowledges 
support from the Research Corporation though a Cottrell Scholarship. 
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